Coupling mechanism between microscopic two-level system and superconducting qubits 
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We propose a scheme to clarify the coupling nature between superconducting Josephson qubits and micro- 
scopic two-level systems. Although dominant interests of studying two-level systems were in phase qubits 
previously, we find that the sensitivity of the generally used spectral method in phase qubits is not sufficient 
to evaluate the exact form of the coupling. On the contrary, our numerical calculation shows that the coupling 
strength changes remarkably with the flux bias for a flux qubit, providing a useful tool to investigate the coupling 
mechanism between the two-level systems and qubits. 
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Recent progress on superconducting qubits suggested that 
superconducting Josephson circuits are promising candidates 
for practical quantum computating (jJ-lSI]. However, extensive 
works are needed to understand the decoherence mechanism 
hence increase the decoherence time of these macroscopic 
quantum systems. For instance, microscopic defects are ubiq- 
uitous in solid state devices. Each of these defects may be- 
havior empirically as a quantum two-level-system (TLS) with 
characteristic frequency ranging about several gigahertz. The 
anticrossings resulting from the resonance between the TLSs 
and qubit were observed in the spectra of phase qubit [6, 7], 
flux qubit [8], and charge qubit BSD. Experiments have shown 
that TLSs not only shorten the decoherence time II 1 Oil . but also 
reduce the visibility of Rabi oscillation of the qubit, limiting 
the fedelity of the quantum gate. Moreover, an ensemble of 
TLSs with various characteristic frequencies may produce low 
frequency 1 // noise ifTTIl . Therefore, it is imperative to under- 
stand the microscopic coupling mechanism between TLS and 
superconducting qubits. 

Unfortunately, it is nearly unattainable to directly probe a 
single TLS's microscopic mechanism due to its microscopic 
nature. Nevertheless, utilizing qubit as a detector of TLS 
supp lied an alternative method to extract useful information 
[l3h . Based on the experimental results, three cou- 
pling models between TLS and superconducting qubit were 
suggested: critical current fluctuator [6], electric dipole lfl2h . 
and flux fluctuator [ 14] . Great effort has been put to determine 
which is the exact coupling mechanism Ill^[l4l4l7ll . Recently, 
it is suggested that one may clarify the coupling mechanism 
by investigating the longitudinal component of the interac- 
tion Hamiltonian d^O]. Lupa§cu et al. have measured the 
flux qubit at symmetric double- well potential and only found 
transverse term [11711 . Furthermore, Cole et al. have theoret- 
ically shown that for different models, the coupling form be- 
tween qubit and TLS changes remarkably. For the electric 
dipole model the coupling type is totally transverse, while for 
the other two models, besides the transverse, a weak longi- 
tudinal coupling exists [14]. However, their experiments in 
a flux-biased phase qubit is unable to determine the coupling 
type, because the observed longitudinal coupling strength is 
too small to discriminate it from the experimental uncertainty. 
Therefore, two questions come out from their works. One is 
whether we have to consider both transverse and longitudi- 
nal terms simultaneously for all superconducting qubits cou- 



pled with TLS. The other is whether the magnitude of the two 
terms depends on the coupling model thereby we can deter- 
mine the exact coupling model by measuring the two terms. 
In this work, we have analyzed the interaction Hamiltonian of 
the three models. We found that the transverse and longitudi- 
nal terms vary for different types of superconducting qubits. 
For phase qubits the longitudinal coupling is always small 
thereby it is not a good system to probe the coupling form. 
However, the longitudinal term for flux qubits is comparable 
to the transverse term. In addition, for different coupling mod- 
els the longitudinal term exhibits distinct flux bias dependen- 
cies, supplying a hopeful scheme to clarify the microscopic 
model of the TLS. 

We start from a short review of the three models which are 
used to describe the microscopic nature of the coupling be- 
tween TLS and qubits. Although they are also valid for flux 
qubits, we at first discuss them in a flux biased phase qubit for 
simplicity. The first model is critical current fluctuator. In this 
model, the microscopic TLS was assumed as a critical current 
fluctuator whose two states respond to two different critical 
currents of the Josephson junction [6]. The fluctuator could be 
considered as a ion moving between the left and right well in 
a double well potential. If fis is the energy difference between 
the two position states, A is the tunneling matrix element, the 
interaction Hamiltonian can be written as II 1 411 

Hi = Vi cos 0(cos 6cr j + sin Ocr^) 

with a X j Z being Pauli operators of TLS in its eigenenergy 
space. V[ = -3^, where 6Io represents the difference of 
the critical currents for the ion populating the right and left 
well, 0o is flux quantization. is phase difference across the 
Josephson junction. 6 denotes the relative orientation of the 
TLS's configuration basis and eigenbasis, tan# = g/A. Al- 
though in many literatures it is assumed = @, Il8ll . we 
would like to keeping the above formula which is more accu- 
rate and general. 

The second model is the interaction of an electrical dipole 
with field [12]. TLS has been modeled as an electric charge 
distribution that changed between the two states. In Ref. 11211 . 
each TLS was assumed as an electron of charge e at position 
R or L separated by a distance d, resulting a dipole moment 
ji = ed. The interaction Hamiltonian is given by 

Hi = ~/3 - E = v q q(cos Ocr x T + sin Ocr^) cos r] 
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where v q = x is the thickness of the junction, q repre- 
sent the number of Cooper-pairs tunneled across the junction. 
7/ denote the relative angle between electric dipole ~jl and the 

electric field E. The defination of 6 is terminologically the 
same as that in the last model. However, s and A may de- 
pend on different physical variables. It is worth to note that 
the coupling term in this model is purely transverse no matter 
whether the double well of the TLS is symmetric ffT^L IT4L ITrh . 
This is a useful characteristic which offers a possibility to dis- 
criminate this coupling mechanism from others. 

Another reasonable model for TLS is magnetic flux fluc- 
tuator. It has been found nearly three decades ago that the 
flux embraced by a rf- SQUID fluctuates with frequency lying 
in low frequency regime. The flux fluctuation is responsible 
for the dephasing of superconducting qubits. Recently, sev- 
eral experiments were carried out to probe the flux noise and 
found that the spectral density of the low frequency flux noise 
with the form l/f jlSWAM . Moreover, Shnirman et al. have 
shown that a collective of TLS's with a natural distribution ac- 
counts for both high frequency and low frequency 1 // noise 
II lh . Therefore, a new explanation for TLS in superconduct- 
ing phase qubit was proposed, suggesting that the states of the 
TLS may modulate the magnetic flux (p e threading the super- 
conducting loop full . The resulted coupling is on the variable 
for the phase qubit. The interaction Hamiltonian can be 
written as 
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FIG. 1 : (Color online) Schematic diagram of the spectrum of a phase 
qubit resonantly coupled to a TLS. The green dashed line and red thin 
solid line represent the resonant peaks of the two-photon transition 
from state 1 to state 4 of the composite system without or with lon- 
gitudinal coupling, respectively. Without longitudinal coupling, the 
two-photon spectrum locates at the exact middle of 1 — > 2 and 1 — > 3 
one-photon transition lines. Note that with longitudinal coupling the 
two-photon spectrum can move up or down relative to the middle 
position, corresponding to positive or negetive sign of the longitudi- 
nal coupling. Here we chose positive sign as an example. Inset: the 
energy levels of the resonant qubit-TLS system. 



Hi = z^0(cos#<x£ + sin#crp 

where = -^(^) 2 . d(f>e is the difference of e between 
the two configuration states of TLS. 6 is defined the same as 
before, with s and A relying on different physical variables. 

The experimental investigation on the coupling mecha- 
nism is not very convincing. Recently, Cole et al. proposed 
a scheme to probe the coupling mechanism in a phase qubit 
(3. Therein, whether the longitudinal coupling exists is a 
crucial clue to decide which model is true. The longitudi- 
nal coupling between resonant qubit-TLS leads to asymme- 
try of the two-photon transition spectrum relative to the one- 
photon transitions spectra: oj\^4 ± co\^i + £^1^3 (1-4 denote 
the eigenstates of the coupled TLS-qubit system), as shown in 
Fig. 1. Therefore, one could experimentally resolve the cou- 
pling type of qubit-TLS system via spectral analysis. How- 
ever, their experiment can not confirm the existence of the 
longitudinal coupling because its value is comparable to the 
measurement uncertainty. Hence, they could not reach a con- 
clusion on the correct model. In our opinion, the key reason 
for the ambiguity is that the longitudinal is much smaller than 
the transverse coupling in all three models. To demonstrate 
this, we unify the interaction Hamiltonian 

Hi = v k O(cos Ocr^ + sin#<x^) 

where k = i, 0, q and 6 = cos 0, 0, q. Ignore the mixed terms 
such as cr z q <j x T , which have no effect on the spectrum of the sys- 
tem, we can write the interaction Hamiltonian in the eigenen- 
ergy basis 

Hi = v k (o k x cos Ocr x q crj + o\ sin Ocr^cr^) 



where the factors o\ z are given by 



ok= Kl|O10) + (0|O|l)l 



o k z = 



Kl|O|l>-<0|O|0>| 



For the electric dipole model, q has no diagonal elements, 
0, so there is no longitudinal coupling in this case. Turn- 
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ing to the other two models, we numerically calculated d 
as functions of flux bias using the parameters in Ref. [I14ll . 

shown in Fig. 2. It is found that < Moreover, for 
(p e > 0.6, the qubit can not work due to the large tunneling 
rate of the excite state; for cp e < 0.57, of is at least one order 
of magnitude smaller than <9^.The much smaller longitudinal 
coupling factor relative to the transverse one may be the main 
reason for which one can not verify the existence of the lon- 
gitudinal interaction between qubit and TLS. Even worse, the 
corresponding coupling factors of the two models are roughly 
equal. This is easy to understand if we notice that in phase 
qubit, - I • We substitute with | + <p, where (p is a small 
quantity, 
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I - into the coupling factors expressions, we 



Taking cos 

can obtain: o l x (f x , o\ ^ o 9 z . Therefore, for phase qubits, the 
longitudinal coupling is not sensitive to the coupling mecha- 
nism. Then, it is difficult to clarify the coupling nature be- 
tween TLS and phase qubit . Although one may argue that the 
TLS parameter 6 has a crucial effect on the longitudinal cou- 
pling magnitude, unfortunately, till now, people are unable to 
control the angle 6 due to the poor knowledge of TLS. 
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FIG. 2: (Coler Online) Coupling factors of qubit- TLS interaction 
in phase qubitb for the critical current fluctuator and flux fluctuator 
models. The parameters of the phase qubit are: C = 850/F, L = 
720/?//, I = 9S4nA. The green dash-dotted and cyan dashed line 
denote the transverse and longitudinal coupling factor o l x ,o[ of the 
critical current fluctuator model respectively. The red solid and pur- 
ple dotted line denote o%, o\ of the flux fluctuator model respectively. 



the upper two levels, showing a very good nonlinearity which 
enables the spectroscopic experiment will not involve the third 
level of the qubit. 

In the three-junction flux qubit, each junction has the possi- 
bility of containing TLS. Even though, we can prove numeri- 
cally that the location of the TLS in different junctions would 
not affect our results qualitatively. Therefore, we consider that 
the TLS is in the small junction without losing generality. For 
electric dipole model, the interaction term is 

Hi = jl-E 

de ~ 

= — V(cos 6crl + sin 6o- z T ) cos tj 
x 1 

_ cte <^o 20 m ( cos Qg-x + sm Qcr 1 ) CO s rj 

X Z7T 1 

Using Heisenberg equation m = ^[0 m , H q ],we obtain 
tia) q d(0\$ m \l) 



Hr 



■ COS 7] COS 6 • CTqCTj 



(2) 



where hoj q is the eigenenergy of the flux qubit. Obviously, 
similar to that in phase qubit the longitudinal coupling is zero. 

For the other models, following the same procedures used 
in the phase qubit, we can straightforwardly write out the in- 
teraction Hamiltonians. 
Critical current model: 



Instead of phase qubit, we find that flux qubit is a possi- 
ble system to reveal the coupling nature of TLS and qubits. 
We start from a flux qubit which consists of a superconduct- 
ing loop interrupted by three Josephson junctions [122ll . Two 
junctions are the same and the other one is a times smaller 
than them. If we assume that the large junctions have a criti- 
cal current Iq and a capacitance C, then the critical current and 
capacitance of the small junction are alo and aC, respectively. 
The Hamiltonian of the circuit is 1221 



H T 



2n 
2n~ 



cos(27r/ + 20 m )(cos 6cr j + sin Ocr^) 
(p\ cos 6cr x q cr x T + o\ sin 60^0*.) (3) 



o, = 



|(1| cos(2tt/ + 20 m )|O> + (0| cos(2tt/ + 20 m )|l>| 



|<1| cos(2;r/ + 20 m )|l> - <0| cos(2;r/ + 20 m )|O>| 



4E 

H q = AE c h\ - Ejcos 0i + 4E c hl - Ej cos 02 + — -n\ 

-aEj cos 03 + Ejil + a) (1) 



where E c = jg, Ej = 7o0o/27r. fa (i = 1,2, 3) is the phase 
difference across each junction, and its conjugate variable hi 
is the number of Cooper pair through each junction. If the 
external magnetic flux <p ext = f<po, using the flux quantization 
condition, we get 03 = 2nf + 0i - 02. Transforming the co- 
ordinates 0i, 02 to the sum and the difference of the phases, 



+ 02)/2, <p n 



■ 02)/2, H q is reduced to 



Magnetic flux fluctuator model 

Hi = 2naE jS(p e sin(27r/ + 20 m )(cos 6cr^ + sin 6cr z T ) 
= 27iaEjS(p e (oi cos Q&ffi + of sin Qo\o\) (4) 



Or = 



Oy = 



|<1| sin(2;r/ + 20 m )|O) + <0| sin(2;r/ + 20 m )|l)| 



|<1| sin(2;r/ + 20 m )|l) - <0| sin(2;r/ + 20 m )|O)| 



H q = E p h 2 + E m h 2 m - IE j cos P cos m + Ej(2 + a) 



-aEj cos(2^-/ + 20 m )) 



where E p - 2E C 



- -i— E 



E p /(l+2a), 



We have calculated the lowest three energy levels near / = 0.5 
, with typically chosen parameters Ej/E c = 40, a = 0.68. In 
the region 0.49 < / < 0.51, the energy difference between the 
lowest two levels (qubit) is much smaller than that between 



Now we compare the magnitudes of the transverse and 
longitudinal couplings. As discussed before, currently it is 
impossible to change the orientation of the TLS basis, we fo- 
cus on the factors d^ z . Using the same parameters as above 
(Ej/E c = 40, a = 0.68), we have numerically calculated o l £ z 
as functions of /, shown in Fig. 3 and 4. When / varies 
from 0.5 to 0.51, d£ z exhibit remarkable changes with totally 
different trends. For critical current fluctuator model (Fig. 
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FIG. 3: (Color online) Transverse (solid line) and longitudinal 
(dashed line) coupling factors for critical current fluctuator model. 
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FIG. 4: (Color online) Transverse (solid line) and longitudinal 
(dashed line) coupling factors for flux fluctuator model. 

3), the transverse coupling factor o\ is zero at the degener- 



ate point / = 0.5 while the longitudinal factor reaches the 
maximum. Then, away from / = 0.5, the transverse factor 
increases monotonically while the longitudinal one decreases 
gradually. At / = 0.51, the transverse factor becomes much 
larger than the longitudinal one. In recent experiments 111 711 . 
the splitting resulted from transverse coupling is observed at 
the degeneracy point of a flux qubit. In addition, no longitudi- 
nal coupling is observed at / = 0.5. These behaviors disagree 
with the predictions of the critical current fluctuator model, 
indicating that the critical current fluctuation is not the domi- 
nate mechanism of the qubit-TLS coupling. 

For the flux fluctuator model, the trends are totally converse 
(Fig. 4). The transverse factor reaches maximum at the degen- 
erate point while the longitudinal magnitude vanishes, indicat- 
ing that the coupling is pure transverse at this point. However, 
the electric dipole model predicts similar pure transverse cou- 
pling [see Eq. d2])] . we can not discriminate the flux fluctuator 
and the electric dipole model at the degenerate point. When 
we tune the flux bias away from / = 0.5, the transverse fac- 
tor decreases and the longitudinal one increases gradually. At 
/ = 0.51, the longitudinal factor is larger than 1 .5 times of the 
transverse one. Therefore, the flux fluctuator model contains 
both transverse and longitudinal coupling but in the electric 
dipole model only transverse interaction exists. We can clar- 
ify the microscopic mechanism of TLS by studying the cou- 
pling term of TLS-flux qubit interaction in a flux qubit biased 
away from the degenerate point. In practical, TLSs have been 
observed in three- junction flux qubits biased away from the 
degenerate point [8|], suggesting that this spectral method is 
completely feasible with the current technique. 

In summary, we have calculated the qubit-TLS coupling 
factors of both transverse and longitudinal terms under three 
microscopic models. It is found that for phase qubits the 
longitudinal coupling is difficult to observe because it is al- 
ways much smaller than the transverse coupling. Then, we 
show that in three-junction flux qubit the relative magnitude 
of the transverse and longitudinal coupling factors are largely 
model-dependent and very sensitive to the external flux bias. 
We propose that these features can be used to clarify the mi- 
croscopic model of TLS. 
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